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ABSTRACT 


Field patterns are drawn for a plasma filled waveguide of circular 
cross section immersed in an axial magnetic field. The plasma is assumed 
to be cold and collisionless, and linearized field equations are used. 

Besides the exact analysis the MHD approximation (low frequency) and 
the quasi-static approximation (slow wave) are used. Comparisons are 
made among field patterns determined by the various methods of analysis. 
Among the field patterns at different frequencies and wave numbers 
certain features are pointed out and comparisons are made. These field 
patterns include particle charge density and current density, total 
charge density and current density, the electric field and magnetic 
field. 
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FIELD PATTERNS FOR A PLASMA FILLED WAVEGUIDE 
LI. INTRODUCTION 


The objective of this thesis has been to draw field patterns 
for a specific plasma filled waveguide in order that field patterns 
of plasma filled waveguides in general may better be visualized. At 
the beginning of Chapter 1 several assumptions are set forth. From 
these assumptions the necessary equations for the fields are derived. 
In Chapter II, the equations are organized in a manner used in arriv- 
ing at the field patterns and then specific regions of the ( 3-W ) 
diagram are discussed. The ( @-w ) diagram is divided into the 
magnetohydrodynamic (MHD) region, the cyclotron resonance region, 
the plasma resonance region, the region of cutoff fields (propaga- 
tion constant equals zero), and a region with the propagation 


constant intermediate between cutoff and resonance. 


° 
too 


te 


. 





II. EQUATIONS GOVERNING FIELDS 
In this chapter the equations of the fields are derived. Eight prelim- 
inary assumptions are made concerning the plasma model and the type of 
geometry. The first derivation is the dielectric tensor which comes from 
applying conservation of momentum principles to the particles. Under the 
heading of the determination of the fields the transverse fields are found 
in terms of the axial fields and the dispersion relation is derived. 
Under the heading for the determinantal equation boundary conditions are 
matched for fields having no variation in the angular direction. For this 
case the fields are written out in explicit form. Under the heading of 
cutoff fields a method of solution is given for finding the cutoff frequen- 
cies. A discussion of the MHD-regime is then made and finally the theory 
of the modified quasi-statics approximation is developed. 
2.1 Preliminary Assumptions | 
The field plots of this thesis are for a very restricted situation; 
therefore the development of the equations will be restricted from the 
beginning and further restricted as the development of the equations 
progresses. A more general development which this thesis parallels 
closely is given in (1.) Waves in Anisotropic Plasmas, PP. (18-24), 
CH. 9, and CH. 10. The preliminary assumptions and restrictions are: 
(1) A linearized macroscopic description has been chosen for 
the analysis. This means that every time dependent field 
Quantity that is determined is a small perturbation from 
an unperturbed state. 
(2.) The system, waveguide and plasma is homogeneous in time 
and space. Fields of the form ei 4 t wi00 be wedtined 


and the plasma characteristics do no change with position. 


(39) 


(4. ) 


(5.) 


(6.) 


(7.) 


The geometry will be restricted to cylindrical coordinates 
with the radius of the waveguide a= .0O5mn. 
The conducting surfaces of the waveguide are assumed to be 
perfect electric conductors, i.e., n X EO atre=a, 
The plasma description is further simplified by assuming that 
no collisions occur and that the plasma is temperate. A 
temperate plasma is one which satisfies the inequalities that 
the particle induced velocity be much less than the random 
thermal velocity* which must be much less than the phase 
velocity of the wave. The lowest phase velocity which is 
considered is in the magnetohydrodynamic region where the 
phase velocity is as low as 1.75 X 10°© times the velocity 
of light. For the electrons to have a thermal velocity 
less than one-tenth the phase velocity of the wave in the 
MHD region the temperature of the plasma must be less than 
10+ °K or in terms of electron volts the plasma must be 
less thana 1 eve plasma. All the rest of the patterns 
have a phase velocity greater than one-tenth the phase 
velocity of light. 

*For a Maxwell velocity distribution the most probable 


velocity of a particle is given by up =, 2kI where k 
m 


is Boltzman's constant, m is mass, and T is temperature. 
It is by definition the thermal velocity. 
The axial unperturbed magnetic field has a strength of .1023 
webers / m-. * 
The plasma is made up of electrons and ions, i.e., protons 


which have a mass of 1.836 X 103 times the mass of an electron. 





(8.) The unperturbed density of each type of particles is 

1.808 x 10L7 particles/m>. * 

*The values of the quantities in assumptions(6.) and (8.) were 
chosen as such for two reasons. The first is that a 
normalization of quantities which is used in the derivation 
of equations gives convenient numbers to work with. 
Secondly, the order of magnitude was chosen because they 
are thought to represent typical values for a plasma filled 
waveguide. 

2.2 Determination of Dielectric Tensor 
The equations governing the fields are Maxwell's equations and 


a conservation of momentum equation for each type of particles. 


VX BF -Jwio (1.) 
VX Ht jweok fs (2.) 
jum, ve = e (Ef vy X Bo) (3.) 


where the k subscript stands for the type of particle. 
Since there is uniformity in the z direction, the fields will be 
broken into their transverse (T) and axial (z) components with Bo = Bo iz 


Solving equation (3.) for v gives _ 
|jw- esBo i, xX 


Mug * i e 
) k E 4. 
PT eeBey | ay PY 
Mk 
“ez “ Sk Ez (5.) 
Jw k 


In order to have the units stand out some normalization constants 


will be introduced and some quantities will be defined. 





Wy, =F (6.) 
aes a (ae) 
k 


ati 





Yk = RS oa (8. ) 





GM wre 

u = UJ (9.) 
Wn 

a2 2 

uy We Ft UD (10. ) 


where c is the velocity of light in free space and a is the radius of 
the waveguide. u,, is the normalized cyclotron frequency, Upk is the 
normalized plasma frequency, Up is the effective normalized plasma 
frequency, and u is the normalized frequency. Assumptions (6.), (7.), 
and (8.) gives ug = —3 and Up = 4. 

Substituting the normalization constants into equations (4.) and 
(5.) and then writing the velocity in matrix form gives 


— ju —u 














| ck | ‘ 
oe 
Vk sk Uck _ oa 7 1 Oo -E (158) 
Mm Wry aa wok | ue _ = | 
0 i... lL 
| | - 


In order to determine SJ in terms of EB some definitions are 


introduced. 
Ax = No & Nx (12. ) 
A 7 de + Si (13-) 
JZ = f+ £ (14, ) 
T 
kK ole t oa (15.) 


is determined from equations (11.), (12.), (13.) and (14.). 
From equations (12.), (13.), and (14.) 
GE = * (en ve + ey Wy) (16.) 
Substituting equation (11.) into equation (16.) determines @. The 
@ thus found substituted into equation (15.) determines K the dielectric 


nw 


tensor. Carrying out the above operations give 


9. 


= 








K | x | O 
ay xX 
ie Meee =k (17.) 
K = Kx | K | 0 . 
oN et —- — +— — - 
0 T O | K 
| i 
where 


acs (ue + Uce YUoi) 


=l— : 
Ky (ae) (u2— 2 (Leas) 
. —Jju cor (uaa + US) a 
x = 19. 
(ue— W ce ) (ue aoe) 
e ; oe 2 (20. ) 
, 1- (2) 


Waves in Anisotropic Plasmas has the equations of this section on 
pages 18 - 23. Equation (11.) is equation (2.8) on page 19, equation 
(17.) is equation (2.11) on page 20, and equations (18. - 20.) are 
equations (2.21) on page 23, 

2.3 Determination of the Fields 

Substitution of equations (14.) and (15.) into equation (2.) 
leaves two vector equations with two vector unknowns. Equations (1.) 
and (2.) are rewritten along with two auxillary equations found by taking 
the divergence of the two curl equations. These auxillary equations are 


not independent of the first two, but nevertheless, are useful relations. 


VX E = “jw H (21.) 
VX H = jwéo kK: E (22.) 
WY Yese = 0 - (23. ) 
VW > gh = 0 (24. ) 


Again the procedure is to split the fields into transverse and 
axial components. Since uniformity exists in the z-direction, the axial 


oz 


dependence of the fields can be taken as e > The operators become 


10. 





Vo ~ yt (25.) 


RK. E = K Ep +k, T,X Bp + Ky cee (26.) 


After splitting the fields into their components, the trans- 


q 
i 


versely directed part of equations (21.) and (22.) can be solved 

for the transverse fields in terms of the axial fields. This is 

matrix (27.). Using this matrix and the remaining parts of equations 

(18 - 24) two coupled equations between Ez and H, can be found in 

a straight forward manner. These are equations (36.) and (37.). The 

above procedure is carried out in more detail in Waves in Anisotropic Plasmas, 
pages 134 - 138. Matrix (27.) is the same as matrix (9.21) on nage >. 

and equations (36.) and (37.) are the same as equations (9.30) and 


(9.3L) repsectively on page 136. 


P | | lg E 
AD A ss ei ee 
Ey Teer iemue || OQ Vp 4, 
E wa. ine 14, | | 
ey Sas) tit) Pst 5) ee 
XR, X Hp =~ | 5 | t TX Vek 
ee 
where 
nee (We acces + (u?x,)° 
a (28. ) 
"= al (29.) 
P = -T(r2 + uk, ) 
D (30. ) 
R = jus K, Mo 
D sae 
Eo (sr) 
Q= PuPKy 
D (ean) 


Ss ju (M2 + ueK, No (33.) 
Eo 
O 


T= gPFu Kx (3h) 


Eo (35.) 


2 He: ex) 
Ky = 


D 
U= -juky (r* + ue Kyt 
D 








The two coupled equations are 








2 
Ve E + ¢ B= a HL (36; ) 
Va. as ae: E (ie) 
Ts Z 8 % ‘ 
where 
caecs pe y eK Ku 
tcl ae. (38. 
ule 
dad = j ul Kx 
a°K | Eo (39.) 
f = (re + ue kK, + ore) 1 
eye aé (40.) 
gs eg ul K 
= St fo (42. ) 
a Ky M, 


Equations (36.) and (37.) can be combined to give two 


- uncoupled fourth order equations. The equation for I is 
2 
ivi + (c+£) Um + (cf - ag)| E, = 0 (u2.) 


Plane wave solutions of the form exp (-jp-xT) satisfy equa- 
tions (36.), (37.) and (42.) where P is a complex function of 
frequency only. The application of the identity (43.) to equation 


(42.) gives the dispersion relation for the medium equation (44.). 


eae ett 
= ° = = saoeEl 
Vp & ep” Pe (43.) 


h 
p -(c +f) p* + (cf - dg) = 0 (4b. ) 
In order to satisfy the boundary conditions at the waveguide wall 


two two dimensional Helmholtz equations will be solved using 


equation (42.) to give the relation between the two P valves. 


a> => 





(Va + p35") ae = «(CO (45.) 


E 


2= Egy + Exp (46. ) 


Ho = D4 Ez4 + ho Exo (47. ) 


4 


where h, is determined from equation (36.) 


hy = Cc - PT 
a 
The above equations can be found on pages (154.) and 


(48. ) 


(155.) in Waves in Anistropic Plasmas. 


2.4 Determinartal Equation for Cylindrical Coordinates 


For cylindrical coordinates equation (45.) is 





>°E 33 + 1 0 Ei + 1 % Bei + Dee Eo4 ears 
a = Ss Bee (49.) 


Applying a separation of variables, equation (50.), to 
equation (49.) gives equations (51.) and (52.). m, is the 
separation constant and is not to be confused with the m 
used for mass earlier. The m for mass will not be used in any 


of the following discussions. 


Boi = Aj R(x) 0 (>) (50.) 


a° 
“age + n, ® = O (51.) 


(pyr)* aoR(r ) 
2+ p.r adaR(r ( = On 
A tpyeh Pt See) * heat) -erve 0 G52) 
In order to have the boundary conditions satisfied for all 
d my = m, and for the fields to be periodic with @ m must 
be an integer. Finally, R(r) = d,@®4r) since there is no need 


for singular solutions in the waveguide. Therefore, the solution 


13. 


Sought is 
Z jmp 
Bi = Ay Um (Pur) e (53.) 
For the rest of this section it is assumed that m = 0. 
For m #0 see page 226.in Waves in Anisotropic Plasmas. The 


boundary condition E, = 0 at re=a gives 


Zz 
A, J (Xo) = ~A,J,(X;) (Su. ) 
where 
Xs Ss Py a (55) 


The boundary condition that Ey = O at r=a gives the 
determinantal equation. E@ is determined from matrix (27.) and 


equations (46.), (47.) and (53.). The determinantal equation is 





X, Ah(X1) _ Xe 4) 
-t2. =; (56) 
bl ox" l bo Jo (Xp) 
where bee = = pe + uc K | 
TM ion 
a D 


Rewriting equation (44.) in its explicit form gives 


ps + ~— (58. ) 
; = s (e +=") — | (c - fp 4dg (59.) 
p 


2 


The simultaneous solution of equations (56.), (58.) and (59.), 
for varying values of u determines the ( R-W) diagram. Fig. (1) is 
the part of the diagram that is complete for m= 0. This diagram 
results from the work of Messrs. P. Serafim and P. Chorney. 

Using the following four definitions, the complex functions for 


the — and H fields for m= 0 are given by equations (64.) through 
(69.). 


14, 





7 (60. ) 
F A Jo (x) . _ | ¢ {fy 2 
; (? (Xp) : 0 (% a. TH (61. ) 
_ oF 
G) a " (62) 
Go = a oFo 
or (63. ) 


2 
Be = D> G, [(P-v2k, lr? + u® x,)+ uw? x, x,? 
i —P(r* + Ue x, +f" K\ (64. ) 


é (Pux,)® +(r? + vK)° Kk, - Kir Peuhe, x, as 
5 ae . ° 
a > 2. ,¢  » 2 1s ee 
e > . Rise. [re * + uw K,) +(u fa) 
15d 
E, = E+ &, (66. ) 


Il 
~juk,, (r + ur K)) +(u°K, } | 





2 
(rux,}° (7° + u°K,) K, - K (7*+u°k, ) ae (67. ) 





xX 
| 
Nl 
o|™ 
Oo 
: [po 
QR 
+ 
—I 





es) 
oO 
tl 
[1 
OQ 
pte 
1 


= ju (K, + Ky (r+ weK,) tu? KP? - Ky 7 (68. ) 


2 2 eS 2 
(ro +u K) + (vx) 





eC 2 
(P+ uk) K, = K, x4* | (69.) 


AA 
0 
mm . |\,€0 2 
EE Ss 
0 
- ja PK, 


By inspection of equations (65.) and (66.) give 











e ~ |Z (—<— E, (70. ) 


15. 


A more general statement of equation (70.) can be found from 


matrix (27.). It applies to field patterns for which Vp =z, 


is parallel to Vp Ez, . By definition 


weet, =" 1, (71.) 
Ve Bea = FS (E Ee (Teas) 


Matrix (27.) then gives 
a 


a> = 
> ER Hy . fy x (Q Vp Eyy +S Vop Hei) 


A= (73.) 


2 
ce ofp = ig e (2 ae Bi t+ EV oe Hs) (74.) 
Azi 





Therefore Ha = alg (2 a (75.) 


2.5 Cutoff Fields (TM = 0) 

By definition a cutoff wave is one whose propagation 
constant equals zero. When cutoff occurs the coefficients 
P, Q, T, d, and & become zero andequations (36.) and (37.) 


become uncoupled. 


Therefore 
2 
(Vp +P) Be = 0 (16.) 
2 a 2 
RK = (Pet) = u Ki (Tie 
e 2 
(Vp + B, )H, = 0 (78.) 
2 2 2 2] 
xX = (P,, 4) = ue (K, + K*) (79.) 
| eee es 





BE - Waves; 
fe 








BL = Ae Im (X¢ = re EY) (80. ) 
= on 

H = 0 (81. ) 

Ep = 0 (82. ) 
= ~Jja Eo! — 

ya - [5 ek Vr E, (83.) 


The boundary condition that E, = 0 at r= a requires X, 
to be a root of the J, (Xe) Bessel function. The cutoff frequencies 


are given by 


nee = xe. + use (84. ) 
H - Waves: 
a 
S r jm ul Ww t 
Ho = And, (X=) ed™? ed (Wt) (85.) 
E, = 0 (86. ) 
Hp = 0 (87.) 
= jue — 
Aa —" D ae H + 2, XxVp 4, (88. ) 
o h K | 
For the H-cutoffs the boundary conditions become more 
difficult; therefore m = O and m = 1 will be considered 
separately. For m = O V pp H, is in the a direction. The 
boundary condition E er O at r = a requires 
any. ie) = 0 (89. ) 


Li 





Again the boundary condition requires Xj, to be a root of the 
Jy (x) Bessel functions. The cutoff frequencies are given by 
equation (79.) which expands into equation (90.). Figures (2.) and 


(3.) are plots of equation (79.) and nine of the solution points are 


indicated. 
ue, + a] uty + 85 uso + az = 0 (90. ) 
where a, = i + ee , + 2u5, + he (91. ) 
ao = (up - Uge Ugi) Bs ye (Use ee +5) (92.) 
Su ne Uci (Ue Uy ~ us ) (93.) 


2.51 H - Cutoff (m = 1) 


Vr H, = _%n Ab eJd® oJu( yt ) r 5 Xx aay - z 
a o(, a. i a i,. .(94.) 
bx, 2) 
J, \X, = 
Xe 


The boundary condition that E p= 0 at r= a gives 


X, Jy (X,, ee Kx (95.) 
J (Xp, ) Ky 


The simultaneous solution of equation (79.) which is 
repeated as equation (96.) with equation (95.) gives the cutoff 


frequencies for m # 1 


eo. 2 (ue + KO) (96. ) 





Ky 
Figures (2.) and (4.) for low frequencies and figures (3.) and 
(5.) for high frequencies demonstrate one method for solving for the 
cutoff frequencies. Equation (95.) is plotted on figures (4.) and 


(5.). The left hand side of equation (95.) is plotted versus X, . 
18. 





The right hand side is plotted on the same axis to the same scale as the 
left hand side but versus u . From these figures a function of 
(X, vs. wu) is generated by selecting an ordinant on the graph and 
reading the corresponding X, and u, then repeating with a new choice 
of ordinant. Equation (96.) has been plotted on figures (2.) and (3.). 
Also plotted on these two figures is the (Xp vs u) function generated 
from figures (4.) and (5.). The dashed curves are for imaginary values of 
X,- Therefore, intersections of either solid lines or dashed lines 
simultaneous satisfy equations (95.) and (96.) and thus are solutions 
from the Hecutoffs from m = 1l. 

For finding field solutions in the cutoff region equations (83.) 


and (88.) are written in a more explicit form. Equation (83.) for 


m= Oand m = 1 respectively is 
Ag (Xe, u) = A. eduly*) Xe\ |€o 
(= = (97.) 
Hp = J Ag (Xe, u) J, (Xe = By (98. ) 
. JO J X t —_ ae 
Br j Ae (X,, u) e ° 4 | e “ ip (Fo (%e 4 ~ I (Xe sl) | 
2s —— (99.) 
Xe a Xe = 


The explicit form of equation (88.) for m = O and m = 
respectively is 


he (Hy) = my eC at) (gn) [% (100. 


er = *, (%, 4 4) J, (Xp =i Th = 3 v (101. ) 


Ki 


19. 


E - na a)e 9 AE volte $) -(e-Me) ln Et, sen 
™* o (% =) - t ‘€ J, ( = 8 


where the characteristics of ( —j kK, / K, 7 : given Mf Figures 
(4.) and (5.). 

Waves in Anisotropic Plasmas considers the cutoff fields on 
pages 158 - 162 and page 227. 
2.6 Low Frequency, MHD Regime 

By definition the MHD regime is the limit that occurs when the 
frequency is much less than both plasma frequency and the ion 


cyclotron frequency. In this limit 








KL ——_~- IL - uy” 
Yoe Uci (1037) 
kK -—» -ju u.= (Uce +4c4) (104. ) 
poet “cx 
a eee Up 2 
"7 (105. ) 
U 
(c - f) ——» — ee I= 
(“4 = (106. ) 
eee. 2 ie 
- ae —| (107. ) 
a Ky 


In the limit of zero frequency ( K.——» constant), (ii —»z 0), 


A 
and (K,,——> -o0 ). A comparison of equations (106.) and (107.) shows 
that |< Cae )- >> 4 dg | ; therefore from equations (58.) 
and (59.) P) == c and pp = ff. Only a slight coupling exists 
between equations (36.) and (37.). The MHD approximation considers the two 


equations to be uncoupled: 


20. 





For the H-Wave, equation (37.) when (g = 0), the cutoff frequency 


is much greater than the ion cyclotron frequency, hence it is not in 


the MHD region. 


KE - Wave: 
in? 
2 
Vp £, + Pp E, = O (108. ) 
2 OR fa 2c Ky! 
a (payo S(n= + uo RS (109. ) 


ae 


Equation (109.), the dispersion relation, is plotted in 
figure (6.). In this figure, the slope (8n/u ) approximately 
equals 67, whereas a slope of 57.1 gives a phase velocity equal 


to the Alven velocity. The definitions of phase velocity and 


Alvén velocity respectively are 


: (110. ) 


Va = “Use Got C (111.) 
| “Pp 


In the limit of zero frequency Q and T, equations (32.) and 





(34.), approach zero; and P and U, equations (30.) and (35.), are 


given by 
Eo ae —la (112. ) 
———— ae 
Pr Seeuse, 
ej au kK € (113.) 
U —— p> Se _2 
ia + ue Ky Ad 
giving r “ — 
Bp = a a t\ 114. 
aoe 
X.2 K V oe 


e be zZ 
ee = > ae He tz xVrEz ae) 


ele 








Equations (114.) and (115.) were used in arriving at the field 
plots for the MHD regime rather than using actual values for the 
coefficients. The results seem questionable and will be discussed 
in the following chapter. 

2.7 Modified Quasi - Statics Approximation 

Near the resonances it is possible to have waves with arbitrarily 
Slow phase velocities in the assumed plasma model. By definition the 
phase velocity of a wave is its frequency divided by its propagation 
constant. Quasi-statics is an approximation that takes advantage 
of the slow phase velocities. In this approximation the velocity of 
light is assumed to be infinite as compared to the phase velocity. 
From this it follows that the E field is approximately curl free. 
(This statement becomes more clear when one assumes plane wave 
solutions of the form exp.( —j k.3) for Maxwell's curl equation 


for E. 
Pv 





-jkK XX E = JWu EB (116.) 
k 
~ x E = Y My H ct 7) 
k 1s sort | as 

oO 


For finite H the right hand side of equation (117.) approaches 

zero.) In Waves in Anisotropic Plasmas pages (162-169) Professor 

Bers studies these resonances by taking limits of the exact solution. 
At plasma resonance he shows that the BS field approaches zero. 

For this case quasi-statics is a very good approximation. At cyclotron 
resonance he shows that i O) and By remains finite. For this 


case the dispersion relations differ in the two approximations. The 
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limit of the exact solution gives for the normalized p-number 


2 2 
X ee ee aT 2 
- ee (118. ) 


while that found from quasi-statics is 


xe = pe a (119.) 


: 

Before proceeding with the derivation it should be noted 
that the field plots are for a modified quasi-statics. That is 
the E field as determined by quasi-statics is used to determine 
an H field, hence the name modified quasi-statics. 

Waves in Anisotropic Plasmas has the derivation of the quasi- 
Pesessic field solutions on pages (240 - 248) and for the modified 
quasi-statics on pages (253-254). Also equation (118.) is the 


same as equation (9.219) page (164.) in the book. 


2.7.1 Derivation of Field Equations 
For the quasi-static approximation equations (120.) and 


(121.) are assumed to describe the fields. 


VW & E = @) C015) 


Er ao (121. ) 
fw 


Since the field is cure free a potential can be defined 


E= «VO (122.) 
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Splitting the field into its components as previous proceedures 


followed give the following equations from equations (121.) and 


(122, ) 
Ve d + PO = 0 (123.) 
Poe Ye Ky (124. ) 
Ky 
Ep =-Vz (125.) 
EB = X*o (126. ) 


Figures (6.) and (7.) are a plot of equation (124.). The 
boundary condition that E,= Oat rsz=a requires (X = pa) 
to be a root of the Bessel function. 

The x- field is determined from the im field and Maxwell's 


equations. 
VX H = JWEK. E (127.) 
m= (0 (128. ) 


For equation (127.) a particular solution can be found which 
satisfies the right hand side and a homogeneous solution can be 
found which will be added to the particular solution in order that 
boundary conditions can be satisfied. 

Again the fields are split into their components and solutions 
of the form exp.(—jp-rp) are assumed. A few algebraic operations 


on equations (127.) and (128.) give 


eh, 


-=_ 
== AJ 
: 


- - 


a 





(129.) 
Hz = 13 ED (130. ) 

where Y, = —ju | €o K 
oe at (131) 

Ee if 
y = J -Eo (132. ) 
Y, = sju éo Ook CK, (133. ) 

ie Ao ae 


\ 
The homogeneous part is found from equations (134.) and 
(135.) 


 X Eo °@ (134.) 
oe” (135.) 

Therefore Ap ye = 7 (136. ) 
i enya! 0 - G5) 


and the total H - field is 


H=H, + EB (138.) 
Equation (138.) is complete when boundary conditions are 
specified. Since the boundarys are perfect electric conductors, 
the arbitrary constant in the homogeneous part is found by 
lettering Hy = O at r= a, 
After defining some symbols, the Ez and H fields are 


written in explicit forn. 
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Til. FIELD PLOTS 


The section describing the conventions has the equations which 
describe the particle motions listed and a discussion of the two 
conventions used in generating the field patterns. The two conventions 
are (1.) polarization ellipses at various locations to describe the 
field strength and a vector to describe its direction and (2.) spacing 
of lines to determine the strength of the field and direction of the 
lines to describe the direction of the field. All the rest of the 
sections consider the plots for various regions of the (B-w) diagram. 
Comparisons are made between various plots and the interesting 


characteristics of certain field patterns are pointed out. 


3.1 Equations and Conventions Used 

The objective of this thesis has been to use those equations 
which were developed in Chapter I in order to arrive at patterns 
of the fields for selected points of the (@-W) diagram. Figure (1.) 
has indicated on it those points for which field patterns have been 
drawn. In order to get the equations for the field patterms E and 
H were determined from those equations developed in chapter one in the 
various sections. The rest of the field patterns were calculated by 


using the following equations. From equations (14.) and (15.) 


fees 2h | So. (K-11). £ 
Ee a \\ Do (EH) - (149.) 


Equations (8.), (11.), and (12.) give 





U=Uck | ur=us), | 
2 “Uck an , 5 (150. ) 
je (Se | gee) wee! Oye 
a, = oe u~=-Uck k | 
a 0 — oe = x 
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a = a t 34 (ipa) 


- 
Px DW p Ve dy (152. ) 
p - a. Teg (153.) 


Taking the curl of equation (2.) and equation (153.) give 
ip = Via E¢g E (154.) 
The constants for the plasma are given by 


ae fe = _16 = 4 (155.) 
“Pi 18360. - 


Poa “20; == : - ; (156.) 


Equations (151.) and (153.) were used to check the correctness 
of the answers. Equation (151.) was used to check that sum of the 
current densities as determined from equation (150.) actually 
equalled the current density as found from equation (149.). 

Equation (153.) was a second method for calculating the total 
charge density and was used to check the correctness of the answers 
as found by equation (154). This check is equivalent to equation 
(23.) V-K-E = 0. Finally, since the computations were 
performed by slide rule, no more than three significant figures 

can be expected. 

No patterns were drawn for the separate particles except for 
frequencies below (u = 1). This is because nearly all the current 
is from movement of electrons at high frequencies. Therefore, the 
pattern for the electrons is the same as the total current. The pattern 


for the ions is nearly the same as that for the E- field. 
28. 
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The justification of this statement comes from assuming uw >? he 


and looking at equation (150.). 


Je = py | _. 
1 a : -3 eet Yo4 i x z | (157. ) 


The most noticable case for disagreement between the patterns for 





Eeand Ji would be for Ep to be linearly polarized. This would 
produce elliptical polarization for Jip with the ratio of major 
to minor axis equal to 1.6 X 1079. Therefore, it is concluded 
that the term (uy Ty X E) has negligible affect on 

the field pattern and the current density pattern for the ions is 
approximately the same as that for the zg field. 

For the field patterns two sets of conventions were used to 
give two sets of field patterns. The first convention has ellipses 
of polarization at various locations to show the relative strength 
of the fields. In each ellipse a line is drawn from the center to 
the outside which indicates the strength and direction of the field 
at time t = oO when the constant of amplitude for the fields is 
taken as a positive real number. As time increases the arrow 
indicates the direction in which the vector rotates. 

The three components of any field are drawn to the same scale 
unless it is noted otherwise by (X Number) on the pattern. The 
particle current and the total current are drawn to the same scale. 
Finally, the scale of any particular (E,J, or H ) was chosen to 


give convenient sized ellipses of polarization. 
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Only two of these views are necessary since the third can be 
drawn from the first two. The (r - @) view is taken looking down 
the waveguide axis ( #4 z-direction) at 2= 0; the ( r - 2) view 
is taken looking across the waveguide into the minus ~-direction 
at P= 0. The(.-z) view actually is generated from the preceeding 
views. It corresponds to looking into the minus r-direction at 
@ = O and moving the polarization ellipse up or down a given dis- 
tance r. depending whether the ellipse is at plus r or minus r 
respectively when $ = O. 

The second convention indicates the strength of the field 
by the spacing of lines. This convention has been applied to all 
of the transverse fields for m = 1 and to the radial fields 
for m =-0. The ¢d-directed component is divergence free for 
m = O and its fields are readily visualized from the patterns 
of the first convention. The Z-dependence of the patterns was not 
accomplished due to lack of time but can be sketched easily from the 
patterns of the first convention. For an increase in 2 the vector 


rotates in the direction of the arrow. When Z« T §=the 
a” By 
vector has rotated by 90° in space. 
The second convention illustrates two features. It clearly 
illustrates the divergence of the fields, especially at cutoff. 
Laying the pattern for the E-field or for the gy, field rotated 
by 90° in space over the charge density contours shows that where 
the density is the greatest the number of lines starting or ending 
per unit area is also the greatest. The second feature is that the total 


transverse pattern rotates in the counterclockwise direction for m = l. 
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Before proceeding with the discussions of the individual regions on 
the (B-w) diagram, it should be noted that the assumed model of the plasma 
may break down at the boundary. Inherent in the model is the assumption 
of specular reflection. Specular reflection accounts for the fact that 
particles cannot flow through the walls of the waveguide as the field 


patterns indicate. 


3.2 Low Frequency, MHD Regime 

For this region, four sets of fields were calculated, two by 
the MHD approximation and two by the quasi-static approximation. 
Such good agreement occurred between the two approximations that only 
two sets of patterns were drawn, one for m = O and one for m = l. 
The lanewed variation between the two approximations was a factor of two 
for Bry gp» 2nd p- 

For this region and exact analysis was not possible because of 
failure to find a solution point on the (A-wW) diagram. A comparison 
of the coefficients which were used with those from an exact analysis 
using the same frequency and propagation constant will indicate the 
degree of the error. For an assumed U = 17 X 167 and /3, 2 1.2 X 
107° one of the p numbers, equation (58.), is very close to c, equation 
(38.), but for the exact analysis it is not necessarily a root of the 


Bessel function since there are two p numbers. Coefficients which would 


have been used for the approximate analysis to find the transverse fields 


are 
P= j 12.3 (158.) 
6 2.5080 (159. ) 
iH =) GEO (160.) 
Me 1.53 (161.) 
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while the coefficients from an exact analysis are 


pes j 12.9 (162. ) 
Q = -2.6 (163.) 
, - Gre (164. ) 
Us j1.63 (ike ) 


From a Comparison of the coefficients and knowing that there 
are two p-numbers each of which are different from the approximations 
raises an element of doubt as to the value of the approximations. 

Assuming the approximations are a valid representation of the 
fields, it can then be concluded that for both m #= O and m= il 
the fields lie in the transverse plane except for J and J.. The 
strong axial electron current is a result of (Ky —co) as (W—- 0). 
Another interesting feature of this pattern is that the ions move in the 
transverse plane in such a fashion that the divergence of the total 
current is negligible compared to the divergence of either particle .~- 
current, i.e., the particle charge density is approximately 3 X 103 times 


greater than the total charge density. 


3.3 Electron Cyclotron Resonance 

Three sets of field patterns, two for m = 1 and one for 
m = OQ, were drawn for this region from the quasi-statics approxima- 
tion. Here larger p-numbers than those used for other regions were 
assumed for two reasons. The first one was to develop some patterns 
with more variations in the radial direction. The second and most 
important reason. is that the dispersion relation as determined by quasi- 


statics approaches the true dispersion relation only for large p numbers. 
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This can be seen by a comparison of equations (118.) and (119.). The 
model for the quasi-statics approximation assumed the velocity of light 
to be infinite compared to the phase velocity of the wave, in this case 
the phase velocity in the transverse direction, hence large p numbers. 

Two points seem to be indicated from these patterns. The first is 
that some of the lines for Jp and Ht enclose upon themselves. The 
Hp for other regions is completely encircling, cutoff and MHD, or no 
encircling lines as plasma resonance. The second point is that a 
comparison of the fields at electron cyclotron resonance with those 
at plasma resonance shows that the H-field for electron cyclotron 
resonance is about an order of magnitude greater than the H - field 
at plasma resonance when the two fields have similar propagation 
constants and similar values for E,. 

3.4 Plasma Resonance 

This is an interesting region because several comparisons can 
be made. AS a first choice the exact analysis and the quasi- 
static analysis at u = 4.08 will be made. Here it can be seen 
that there is very good agreement between the two patterns and field 
solutions. The difference which occurs is that the E - field becomes 
slightly elliptically polarized in the transverse plane for the exact 
analysis and the strength of the Z-directed field as found by the 
exact method is about twice that from the quasi-static method. 

The second interesting comparison is that at similar propagation 
constants the plasma resonance cames closer to being an E-wave than 
Riek von resonance as won discussed in the preceeding section. Also 
E, for plasma resonance is about twice E, at cyclotron resonance when. 


the two have similar magnitudes for Ep 
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The third interesting comparison is to compare the patterns of those 
fields which have solution points that lie on the curve on the (6-w) 
diagram that originates at = 3.94. See Figure (1.). As the propagation 
constant increases the 2 component of the E and J fields grow while the 
H - field decreases. The transverse A - field grows from nothing 
at cutoff to strong radial then to almost circular while the E - field 
starts from almost circular and decreases to almost radial. The Jp field 
is very similar to the Hn field except that it has a strong radial field 
at cutoff. 

The fourth interesting point is that J, and H have @ polarization 
in the transverse plane which is opposite to that of the E-field. 

The E vector rotates in the counterclockwise direction while J and H 


rotate in the clockwise direction. 


3.5 Cutoff Melds (M= 0) 


3-5-1 E ~ Waves 
All that seems apparent for these cutoffs is that the EH field is 
the only one with transverse components and the charge density is zero. 
They are similar to the cutoffs for Zz - waves in an empty waveguide. 
They differ in that the cutoff .frequency is higher than that for the 
empty waveguide and that the empty waveguide has no particles, hence 
no current density. For the empty waveguide (Ueco= Xe ) of equation 


(8.0) and is a root of the Bessel function. 


34. 





Belee H - Waves 
~ 


mm = O3 
Three field patterns were drawn form = 0. For the cutoff 
at u = .0430 the ion and electron charge density nearly cancel 
each other. From the patterns of Aer and Jip the particle paths 


can be drawn for low excitation, i.e., localized motion of particles. 


(a.) The small ellipse is the path of the ion. 
(b.) The large ellipse is the path of the electron. 

For u = 3.94 the interesting feature is that a nearly 
circularly polarized Ep - field produces a strong radially polarized 
J - field. The reason for this is that UE, approximately equals u Ey 
in equation (150. ) hence (Jy = 0). 

The cutoff frequency u = 6.34 has nothing that particularly 
stands out. A comparison of the three patterns shows that as frequency 
is increased, the E-field goes from nearly radial, to almost circular 
with the stronger field being r-directed and then to almost circular 
with the strongest field being @ directed. For m= O the H wave 
at cutoff in an empty waveguide has the E field entirely in the 
@ direction; therefore none of the patterns are like those of the 
empty waveguide. The higher the frequency the more nearly the modes 


look like empty waveguide modes. 
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For m *® 1] there are four sets of field patterns. A comparison 
of the patterns for U = .04326 and m= lwith Us .0430 and m=0 
shows that for both sets the charge density of the electrons tends to cancel 
that of the ions, that is, Pe = ~ 3.75p and e = 4.75 p. The 


particle paths for m = 1 look like 
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(a.) The smaller ellipse is the path of the ion. 


(b.) The large ellipse is the path of the electron. 


In order to visualize the way the electrons tend to cancel 

the charge density of the ions the second convention is the better. 
A comparison of U = 3.37 for m 2 1 with uw = 3.94 

for m = O shows that both have nearly circular polarization 
ellipses across the cross section of the waveguide, and both have 
the same direction of polarization. 

The last two sets of patterns are for imaginary p numbers. 
The fields for U = 5.07 show how large imaginary p numbers produce 
very strong fields near r = a. This is because a modified Bessel 
function grows rapidly as its argument grows. The fields for 


uoeo= 5.77 are for p numbers that are practically equal to zero. 
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3.6 Patterns Between Cutoff and Resonance 


Two sets of patterns have been categorized under this heading, 
one for u= 4.773 and one for u2 1.00. For the patterns at 
u = 4,773 one of the normalized p-numbers is 8.65 and the other 
is 2.40. The large p-—number causes three variations in the field 
pattern as r goes from zero to a. For this pattern it is neces- 
sary to use the patterns from the second convention because the large 
number of variations renders the first set of patterns practically 


useless. 


The special feature the patterns at U = 1 have is that this 
is the only set which was drawn for complex p-—numbers. For it 
a new complex constant of amplitude was chosen in order that the 
fields would be in either the (r or @ )-directions at time 
t = OQ. Like the cutoff fields with imaginary p - numbers the 


fields are Strongest .near the outside of the waveguide. 
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MHD : 
E-Wave (1.) pp. 169-172 
m=0 
u=17x 107° 
8.. = 1.141 x 10° 
p = 2.05/a 
A(z,t) = aed 


a 
n 


(17x107°w jtel. 1jix107* 2) 


Complex Fields: 


E = A(z,t) (-j795 J,(pr)¥, + Jj(pr)t,) 

H = A(z,t) [-J103.8 J, (pr)fy) 

J = A(z,t) [23.7 J, (pri, + j2.38 Jy (pr Ey 
~j4.99x107 Jo(pr )t,] 

J, = A(z,t) [4+1.273x10°° J, (pr )f_,-j2.25x10° 5 (pr yt, 
-jh.99x10°7 J, (pr)f,) 

J, = A(z,t) (23.7 J, (prdf,, + 42. 27x10° 3, (pr )ty 


~j2.72x107+ 3, (pr), ] 
p = A(z,t) [-53.38x107" J, (pr)] 
A(z,t) [-j1.118x107- J.(pr)} 


i 


Po 
A(z,t) [+j1.-118x107- J, (pr)) 


Ps 





winG < 


Quasi-Static: (1.) pp. 21-246, 253-25) 


m=0 

u = 17x107> 

B= 5 87x1073, = .117h 

p = 2.405/a 

tht) 4,04 (27x10~ Wnt=5. 87x10732/a) 


Complex Fields: 


Note? 


E = A(x,t) [-jhl0 J, (pr)f, + J, (pr), ] 

H = A(z,t) [(-6.4hx10~> J, (pr )-3. 35x10"? £7, 
-J103.9 J, (pr)ty +5 (2.64x107° J (pr) 
“1.142x107°)T.] 

J = A(z,t) [12.20 J, (pr if, + j1.230 Jy (pr) ty 
-j4.99x107 J, (pr )t,} 

J, = A(z,t) [6.56x107> J, (pr)¥, = j1.159x10° 3, (pr), 
~jh.99x10"? J (pr )T,) 

J, = A(z,t) [12.20 J, (pr)f,, + j1.171x10° 3 (pr hy 


=4j2.272x1074 J, (pr )t,} 
0 = A(z,t) [=jl. 746x107! Jo(pr)] 
Pg = A(z, t) [-J5. 75x10" 5, (pr)}] 
p, = Alz,t) [4)5. 75x10 J, (pr) 


Field Patterns Drawn For MHD Only 
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MHD : 
E-Wave (1.) pp. 169-172 
m= 1] 
u = 117x107? 
6 = 1.353x107° 


p = 3.832/a : 
A(b,z,t) = A, J (17x10 Wr tel. 353x100" “2/a) 44 


Complex Fields: 




















—- _ ] Jy (pr) 
E = A($,z,t) [j5.93x10 se > ae 
- a. Ty + Jy (prt) 
) 
H = A(d,z,t) {65.1 —— “t + j65. 1(J, (pr) ce = T4] 
(pr) 
a Jy 
J = A(>,z,t) [-17.63 J,(pr) + 15. i: x ae 
ae ) 
-J(1.779 Jj(pr) + 15.85 zy 
~j.99x107 J, (pr )T, 1 
7 
T, = A(b,2,t) [+(9.49x10 7 - ~ 1.672x10° it 
_ 
+41.672x10° (I,(pr) - Ey 
=j4.99x107 J, (prt) 
J, (pr) 
J; = A(b,z,t) [(=17. 62 J, (pr) ” 1. 519x10° as Ee 





J, (pr) 
~j(1.697x10* Jj(pr) = 1.519x10° Wt, 
~}2. 72x1074 J, (pr 35) 


p = A(b,2,t) [~Jjh.02x107" 3, (pr)] 
p, = Aldyz,t) [-J1.32hx107> J, (pr)] 
p, = Alpyz,t) [+)1.325x107? J, (pr)] 





-54. 


Quasi-Static: (1.) pp.21-2h6, 253-25) 


m=1 

u = 17x107> 

2. = 9.36x1072 8= .1872 
p = 3.83/a 


A(b,2,t) = A, 43 (172x107? Wy t-9. 36x10 “32 /a) _i? 


Complex Fields; 




















= .. ae (pr 9} J, (pr) 
EB = A($,2z,t) (JhO9(J,(pr) - f,, - hoo SF, 
Ae oe ] 
( 
H = A($,2,t) [5.1 - . t. + j65,1(J (pr) - 2 oh 
+31. 652x107° i, (pr fo] 
h : J, (pr) 
J = A($,2,¢) [( 912.19 Ij(pr) + 10.95 = ja. 
-j (1.228 Jy (pr) + 10.95 a, 
-j4.99x10° Jy (pr fo] 
J, (pr) 
J, = A(b,z,t) [-(6. Shx107 o(PT “ er} 156x10° wre 
+41.156x10° (J, (pr) om If, | 
=j4.99x107 J, (pr)? 
(pr ) 
Jy = A(o,z,t) ((-12.19 J, (pr) = 1. 04.8x10° a a i 
( 
+3 (<1. 169x10° J, (pr) + 1.048x10° ——— ~ Ee tegen op 


mo 72x1074 Jy (pr ae) 
p = A(byz,t) [-j2.77x107" 3, (pr)] 
A(>,z,+t) [-39.16x1074 J, (pr )] 
A(>,z,t) 1 39.16x1074 J, (pr )] 


Pe 


Ps 


Note: Field Patterns Drawn For MHD Only 
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wt2o 
P= e, 
Pez 3.3x108 P,, 
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ELECTRON CYCLOTRON RESONANCE: 
QuasieStatic Analysis 


m=0 
u = 2. 
R= 10.08 ; G= 201.6 
a 9-62 
p = ; 
A+). syed (2 Wn ted0.08 Tr) 


Complex Fields; 
A(z,t) [-3.547 J, (pr)f, + J, (@r)t,) 
= 1073xa(z,t) [ (1-642 J, (pr) + 1.937x107 42, ( ar) )2, 
_ 2§2. 05 Jy (pr Ey - §(.898 Jo (pr) 
“ 1.937104 Tenia 
= ACz,t) (. 344 Jy (pr )£, - j.430 J, (pr fy 
-J.354 J, (pr )¥,) 
= A(z,t) [-j2.23x10°7 J, (pr)] 


of | Ol 
tl 


a | 
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D 
I 








G2 
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ELECTRON CYCLOTRON RESONANCE: 


Quasi-Static Analysis 


m=1 


u = 2.) 


Bn 


= 7.01 ; S= 10.2 
wy =e 


Z 
A(>,z,t) = pita rm ®) 56 


Complex Fields: 


E 


ml 


cf 


A($,25t) [§-547(5, (pr) - 


a= ) J, (pr) 
uf, - -5h7 5 








1 
+ J, (pr )t_} 


«3 6 J, (pr) 
LO A($,z,t) { (2.3 J, (pr) + 5.31 <r 


3 
-6.hyx1077 (I. (6r) = an 2 


a ay Ean oi 


-j(1.292 J, (pr) = camel I, (@r))¥,) 
a 





+§(2.95 Jj(pr) = 5.31 





ACb,z,t) [(=-3hh Jy(pr) + . 7H i. 


300), 





+§(.430 Jo(pr) = . 7h 4, 


-§.1937 J, (pr)f,} 


A($,z,t) [-j1.611x10~7 J, (pr)] 
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-—~ 
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ELECTRON CYCLOTRON RESONANCE: 
Quasi-Static Analysis 
m=1 p= T.Qa6 
u = 2.18 
A, = 9-93 = 198.6 
Abo2st) = aod (2-16nt-9.93 =) 44 


Complex Fields; 














a J, (pr) J, (pr) 
E = A($,2,t) [5-707 J,(pr) - “S—It, - .707 AS fy 
+ J, (pr)i,] 
: ) 
H = 1077xA($,2,t) [(-1.425 5, (pr) + 3.38 — a 
I, (Gr) 
4 1. 794x107 #(I, (Br) = i - i. 
3, (pr) <a 7 F 
+5 (1-962 J, (pr) = 3.38 —s-— + 1. 79hx10 s 
~3(1.003 J, (pr) + 179x107 3, ( Br) IE] 
(pr ) 
J = A($,z,t) ((-.308 Jj(pr) + . 732 a. ——)t,, 
= ) 
+3 (42h Ij(pr) = .732 aL iz, 
~§.390 J, (pr)T,] 
p = A(p,z,t) [-J2.63x10°7 J (pr)) 








beer exned 


-o=*) 








le 














14 





ee 








iio 


PLASMA RESONANCE: (1.) pp. 16216 
Exact Analysis 


m=O 
u = 4.08 
a. = 10.53 
Py = 248 
= j21€ = JQ, 


n'y, ¥) = aod (he 08 Wyt-10.53 g) 


Comlex Fields: 
E = A(z,t) [-j2.09x107" (J, (pyr) = 6.35x107° 1 (Qn), 
+ 3.9hx107° (J, (pyr) = 2.h9x107 +1 (ar) )%, 
+ (Jg(Pyr) + Lel27x10~ I, (Qo) Fy] 


TI 
i 


A n,t) (eerexie"4 (3 (Pr) = 2. imo" 1, (Q5r))f, 
+j1.722x10 Aes. 1 (Pi?) + 3, 25x107° I 1 (@r ))Ey 
+j6. 25x1072 (J o(Py") - 9,.89x10 a4 I 9 (Oar ))4, J 


Cf 
i 


= A(z,t) [-8.15x10° “(34 (pyr) = a I 1 (Qr))T, 
+}5.17x10°°(J, (pyr) + 3.9x107! 1 (gr ))Ty 
~j2.08x107* (J, (p,r) + 1.427x107> I oar } 


= A(z,t) [-Jx10"7 (1.948 Jj(pyr) + 3. 76x107° 1, (Qor))] 


OD 
! 





f 














(7 
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PLASMA RESONANCE: 
Quasi=Static Analysis 
m = 0 
u = 4.08 
8, = 8.23 3 (3= 164.6 


(4.08 w te8.23 4) 
A(z,t) = ajo! + i a 


Complex Fields; 


E = A(z,t) [-j2.92x107> J, (pr)f,, + J, (pr)t,) 
H = A(z,t) (5. éox1074(5, (pr) - 1. 045x107 1, (AP))L, 
+jh.20x107+ J, (or Jy 
+56.92x107°(3 (pr) = 8.61x1073 1,( Br))f,} 
J = A(x,t) [-1.322x107* 5, (pr)f, + §9. 72x10 J, (pr )f4 


=j.208 Jy (pr ii 


p = A(z,t) [-j1.58x107” J, (pr)] 
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PLASMA RESONANCE =~ FAST WAVE: 
Exact Analysis 


m = 0 

u = 4.ish 
i 1. 00 
P, = 3.149/a 


Po = j1.87h/a = j a5 
(415k w,t-H 2) 
A(z,t) = A,e 
Complex Fields: 
A(z,t) [j(=.686 J, (pr) + .0927 I, (qor) EZ, 
+ (4393 Jy (pyr) - .0782 Ty (aor) Ly 
+ (Jy(pyr) + .1h6h I5(qor))t,) 


E 


H = 107?xA(z,t) [-1.006 J, (pyr) + .1997 I, (aor))¥, 
+j(.252 J, (pyr) + .0620 Ty (gor) )4, 
+§(.795 Jo(pyr) - 0937 I,(aor))T,] 


J = A(z,t) [-.1717 J) (pyr) + 01545 I, (aor) )z, 
+3 (.0437 Jy (pyr) + 00475 Ty (aor) )¥y 
-j(.20) J (pyr) + .0300 I, (qor))f,] 

p = A(z,t) [-j1077(1.090 J,(pyr) + .0728 I,(a,r))] 








ow 
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PLASMA RESONANCE: 
Quasi-Static Analysis 


m= 1 
u = h.26 
Bar= Ome ; S= je2.), 


j (4-26 W, t=9. 62 =) 


A(>,z,t) = Aje oi? 


Complex Fields: 














J, (pr) J, (pr) 
= 1 iL 
E = A(b,z,t) (§-398(J, (pr) - ~Sr—)2,, - .398 ty 
+ J, (pr )f,) 
es a J, (pr) 
H = 1079xa($,2,t) ((.497 Jg(pr) - .847 —— 
a 
+ 1.129x1074[1_( Br) - 4 ie 
Jy a ~), 1,( 8r) 
@j(.350 J, (pr) = . 84.7 —— e 1,.129x10 4 ar fy 
5(.1977 J, (pr )- 1.129x10 “4 I (62008) 
(pr) 
J = A(bezst) ((.1571 Jg(pr) - 268 LO yt, 
ste : 
-j(.1108 J. (pr) - .268 fy 
-j.1992 J, (prt) 
p = A(p,2,t) [=J1.970x1077 J, (pr)] 








oe 
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CUTOFF FIELDS: 


E-Wave 
m= 0 
u = .67 


A(t) = gest O7 Wnt 


Complex Fields: 

B= A(t) (J,(pr)4,1 

H = A(t) (j-1367 J, (pr)f4] 
J = A(t) [-j.1818 J, (prt, ] 
3 = a(t) [-j.1818 J, (pr )f,] 


A(t) [-39.89x107? J, (pr)%, 1 


om cf 
Ke 
i 








ot 
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CUTOFF FIELDS: 


E-Wave 
m=1 
u = 5.5) 
g g ae 


& 


A(>,t) = a 055-54 Wnt 3i? 


Complex Fields: 
E = A($,t) (J, (pr)f,) 


J 
H# = 1.833x107- A(o,t) [- wer? 


pr 
A(>,t) [-j.153} J, (pr)f,] 


J, (pr) 
pr 








ft, 7 J(5, (pr) - 4] 


a 


p = 0 





a 7 
{or ,4 02-0 Ie ‘ 











wh, 





lOO 


-~101~ 


CUTOFF FIELDS: 


H-Wave 
m= 0 
5 w Babs 


(1-) Ugo, = +0430 A(t) = Ao 0H30Wnt 


co 

B= A(t) J,(pr) (2.72x10° Z, = jh.23 £4] 

EH = A(t) J, (pr )i, 

J = A(t) J, (pr) [-j.619 i + 77.0 £4) 

J, = A(t) Jy(pr) [+J2. 30 E. - 77-0 £4) 

T, = A(t) J, (pr) [-J2.92 Z, + 6. 56x107° ty] 
= A(t) Jy (pr ) (1.83x107"] 

(2.) woo = 3-94 a(t) = aod hwnt 


5 - = cher) J, (pr) [-h.98 Z, = j3.88 £4] 
H = A(t) J, (ort, 
J = A(t) J(pr) [§10}.0 i = helé £43 
p = A(t) Jj(pr) [-3.37x107] 


(3.0) u = 6,3 A(t) = a, o46- 34 Unt 


E= A(t) J,(pr) (310 ZT, = j62h £4) 

H = A(t) J, (pr )i, 

J = A(t) J,(pr) [-j10}.6 a m3 322 t4) 
p = A(t) J(pr) [2.10x107")] 
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